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Abstract 



A unitary quantum gas is a gas of quantum particles with a binary interaction of infinite scattering length and 

^ ' negligible range. It has been produced in recent experiments with gases of fermionic atoms by means of a Feshbach 

>— ' , resonance. Using the Fermi pseudo-potential model for the atomic interaction, we show that the time evolution 

^^ ■ of such a gas in an isotropic three-dimensional time dependent harmonic trap is exactly given by a gauge and 

^T^ ' scaling transform. 

fy^ ' Resume 

f^ ' Nous entendons par "gaz quantique unitaire" une assemblee de particules dont le mouvement est decrit quan- 

tiquement et qui interagissent par un potentiel de longueur de diffusion infinie et de portee negligeable devant leur 

j^ ' distance moyenne et leur longueur d'onde thermique. Un tel gaz a ete produit recemment a I'aide d'une resonance 

de Feshbach dans un gaz d'atomes fermioniques. En modelisant les interactions entre particules par le pseudo- 

potentiel de Fermi, nous montrons que revolution d'un gaz unitaire dans un potentiel de piegeage harmonique 

'■^ ' isotrope tridimensionnel de dependence temporelle quelconque est decrite exactement par la composition d'un 

S , changement d'echelle et d'une transformation de jauge. 

o 
o 
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1. Introduction 

Experiments with quantum gases of spin 1/2 fermionic atoms are currently making rapid progresses. One of 
the most fascinating properties of these fermionic gases is the possibihty to freely tune the sign and the strength 
of the atomic interactions without reducing the lifetime of the sample: the value of the s-wave scattering length 
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a of two particles with opposite spin components can virtually be adjusted from — oo to +00 with a Feshbach 
resonance technique [1] , without inducing any instability of the gas even in the unitary limit a = ±00 [2] . These 
systems are still gases, in the sense that the effective range of the interaction potential is negligible as compared 
to the mean interparticle separation and to the thermal de Broglie wavelength. We shall take advantage of this 
crucial property and model the true interaction potential by the so-called Fermi pseudo-potential [3]. 

Such stability of the strongly interacting Fermi gases opens up fascinating possibilities, e.g. the study of the 
crossover between a Bose-Einstein condensate of dimers (already observed, see [4,5,6,7]) and a ECS condensate 
of pairs, by passing through the strongly interacting regime kF\a\ ^ 1, where kp is the Fermi momentum [8,9]. 
In the unitary limit fcF|a| -^ +00, the thermodynamic properties of the spatially homogeneous gas are universal: 
they depend only on the Fermi energy and on the temperature. At zero temperature, the chemical potential of the 
homogeneous gas is then /i — rjno where r; is a pure number and /lo is the chemical potential of the ideal Fermi 
gas [2,10]. An accurate measurement of rj would provide a crucial test of many-body theories [2,5,7]. 

The standard imaging technique used with quantum gases is to switch off the trapping potential, to let the 
gas expand and to perform a light absorption imaging of the atomic cloud. Such a ballistic expansion acts as a 
magnifying lens: it was used to reveal the vortex lattice in a rotating Bose-Einstein condensate [11], and very 
recently to obtain the value of the universal number rj for the unitary Fermi gas [7] . Clearly the interpretation of 
the time of flight images strongly relies on a theoretical understanding of the time evolution of the gas in a time 
dependent harmonic potential. In the case of a pure Bose-Einstein condensate in the regime kF\a\ <S 1, this was 
achieved starting from the Gross-Pitaevskii equation by a gauge plus scaling transform [12,13]. When the Bose or 
Fermi gas enters the strongly interacting regime kF\a\ > 1, no solution starting from flrst principles is available 
[14] and one relies on the hydrodynamic approximation [15]. 

Here we consider the idealized case of an isotropic and harmonic three-dimensional trapping potential. In the 
limit of an infinite scattering length, we show that the Fermi pseudo-potential has a scaling invariance that 
rigorously allows the use of a gauge plus scaling transform similar to the one of [16] to describe the time evolution 
of the gas due to an arbitrary variation of the trapping frequency. 



2. The model based on the Fermi pseudo-potential 

Consider an assembly of N non relativistic particles, with an arbitrary spin. These particles may be indistin- 
guishable bosons or fermions, or even be distinguishable. All the particles have the same mass m and interact via 
the same binary interaction potential independent of the spin degrees of freedom. The interaction potential is the 
Fermi pseudo-potential with coupling constant g related to the s-wave scattering length a hy g — A-Kh^a/rn. At 
this stage, < |a| < -l-co, we shall take the unitary limit |a| — > -l-oo later. 

Let tpijci, . . . , tn) be the wavefunction of the gas corresponding to a given (but arbitrary) spin conflguration. 
The wavefunction ip then evolves according to the Schrodinger equation: 

'^ r fi2 1 



l<i<j<N 



Here Arj is the three-dimensional Laplacian with respect to the spatial coordinates ri of particle number i, U 
is the trapping potential seen by each particle and each 1/'];^^, the so-called regular part of i/) in ri — r^, is the 
following function oi N ~ 1 vectors of coordinates: 



C;^(K,fc/i,j},Rij)= lim„^(r,,V) (2) 



where r^ is the norm of rij = ri — rj and where the limit and the partial derivative are taken for flxed positions 
Tk of the N — 2 particles other than particles i and j, and for a fixed position of the center of mass of the particles 
i and j, Rij = (n +rj)/2. 

The domain of the Hamiltonian corresponding to the Fermi pseudo-potential is therefore not the Hilbert space 
of the non-interacting gas, but a functional space with specific boundary conditions for the wavefunction ^p. More 
precisely, as we now see, the model amounts to replacing the true interaction potential by contact conditions, 
i.e. by boundary conditions on ^ when the distance Vij between two particles tends to zero, the wavefunction t/) 
otherwise evolving with the interaction free Schrodinger equation [17]. As the wavefunction V does not contain 
any delta singularity, each delta singularity in the interaction term of Eq.(l) has indeed to be compensated by 



a delta singularity in the kinetic energy term. In 3D this implies that ?/> can diverge as l/r^j when Vij ^ 0, a 
divergence which is still square integrable. Two cases can occur: 

- case i): lirur . ^o V' — 0- no delta singularity occurs from the kinetic energy operator, and there is no delta 
singularity from the interaction term as ipl''^ vanishes. For instance, this is the case when particles i and j are 
indistinguishable fermions in the same spin state. 

- case ii): xp has a l/nj singularity: 

V=— + B + 0(r,,) (3) 

Tij 

where A and B are still functions of the rk's with k 7^ i,j, and of Rij. The regular part of ip is then ip'^^^ — B. 
Writing the kinetic energy operator for the pair of particles i,j as Ar; + Ar- — ^Ar;. + 2Ari. and using 
A(l/r) = — 47r(5(r), we find that the total coefficient of S{ri — Vj) in the right hand side of Eq.(l) vanishes 
provided that A + aB = 0. 
A way of summarizing the two cases is then simply to impose the boundary conditions: 

V(ri, . . . ,rN) = A({rk,fc / i, j}, Rij)[r-.' - a'^] + 0(n,) (4) 

the first case corresponding to ^ = and the second one to A 7^ 0. Having ensured that all the delta singularities 
cancel in the Schrodinger equation, we can now restrict it to the manifold where the positions of the particles are 
two by two distinct: 



ihdtip = ^ 



4A,+f/(rO 



rp for Uj / 0, Vi / j. (5) 



Eq.(4) and Eq.(5) constitute the basis of our model. 



3. Scaling transform in the unitary limit 

We specialize the previous section to the case of the unitary quantum gases: the scattering length is now infinite, 
so that we set to zero the 1/a term in the boundary conditions (4), to obtain 

^(r„...,r:.) = ^i«Iilli£MlM+o(,,^). (6) 

Note that the unitary limit does not look in any way singular in our formulation of the model, as the scattering 
length only appears through its inverse in Eq.(4). Also we restrict to the case of an isotropic harmonic potential: 

f/(r) = \mu\ty (7) 

where the oscillation frequency uiit) of a particle in a trap is constant and equal to cjo for f < and has an 
arbitrary time dependence for f > 0. The ballistic expansion of the gas mentioned in the introduction corresponds 
to setting oj to zero for i > 0. 

We assume that the state vector of the gas for f < is a steady state of Schrodinger's equation with the energy 
E. The corresponding wavefunction for the considered spin configuration is ipo, which in particular obeys the 
boundary condition Eq.(6). Our ansatz for the time-dependent wavefunction, inspired from [12,13,16], is then 

V'(ri,...,rN,t)=AA(i)e ^:'=i '^" ^ Vo(ri/A(f), . . . , rN/A(t)) (8) 

where the time dependent normalisation factor N{t) and the time dependent scaling factor \{t) need to be 
determined. This ansatz is the combination of a gauge transformation (see the Gaussian phase factor) and a time 
dependent scaling transform (see the rescaling of the coordinates by \{t) in -i/'o)- 

The first step is to check that the ansatz Eq.(8) is in the right Hilbert space, i.e. that it satisfies the boundary 
conditions Eq.(6). The scaling transform indeed preserves the boundary conditions: it rescales and multiplies the 
function A by A but does not lead to the appearance of a 0(1) term in Eq.(6). Note that for a non-zero value 
of 1/a the conclusion would be different if A 7^ 0. The gauge transform also preserves the boundary conditions: 



for fixed rk's we write r^ + r| — 2Rfj + rfj/2 so that for a fixed Rij, the gauge transform involves only as a 
non-constant factor 

g™r-f^A/4fiA ^ ^ _^ imr%\/4hX + 0{r%). (9) 

As rfj X A/rtj = 0{rij), the boundary conditions are preserved by the gauge transform, a conclusion that extends 
to the non-zero 1/a case. 

What is left is to check that the free particle Schrodinger equation Eq.(5) is satisfied by the ansatz for an 
appropriate choice of \{t) and A/'(t). One calculates the time derivative and the Laplacian of the ansatz. One uses 
the fact that -j/iq is an eigenstate of energy E to express the action of the kinetic energy operator on ipo in terms 
of r^Tpo terms and Eipo- Equating the terms r^ipo on both sides of Eq.(5) leads to 

A(i) = ^-c.^WA(t) (10) 

to be solved with the initial conditions X{t < 0) = 1, X{t < 0) = 0. For a ballistic expansion, one finds A(i) = 
(1 + LOot Y' . Equating the terms proportional to ipo gives Af/jV — — 37VA/(2A) — iE/{hX ) which is readily 
integrated in Af{t) = X~^'^'^{t)exp[—iEL dr/A^(r)ft]. The first factor in the right hand side of this equation 
ensures the conservation of the norm of the wavefunction. 

The ansatz Eq.(8) therefore gives, in the unitary limit, the exact time evolution of the initial state vector for 
an arbitrary time dependence of the isotropic harmonic potential. It rigorously confirms the scaling law that 
one would obtain from the zero temperature hydrodynamic approximation. It also allows to use the symmetry 
considerations developed in [16] : e.g. the limit of linear response to a sudden change in the frequency uj away from 
UJQ gives rise to an undamped oscillation at frequency 2ujo which reveals the existence of A''-body stationary states 
of energy E ± 2^ia;o coupled to Tpo by the excitation operator X^^^j r^. 

We thank L. Pricoupenko for discussions, A. Sinatra, I. Carusotto, J. Dalibard and C. Cohen- Tannoudji for 
comments. LKB is a research unit of Ecole normale superieure and Universite Pierre et Marie, associated to CNRS. 
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